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Parthenogenesis is a natural form of asexual reproduction 
in which offspring are produced from unfertilised eggs 
(Neaves & Baumann 2011, van der Kooi & Schwander 

2015) . This uncommon reproductive strategy has been re¬ 
ported in less than 0.1% of vertebrate species, neverthe¬ 
less including a wide range of taxa (i.e. fishes, amphibians, 
reptiles, birds and mammals), even in wild populations 
(Neaves & Baumann 2011, Booth et al. 2012, van der 
Kooi & Schwander 2015, Dudgeon et al. 2017, Allen et 
al. 2018). Parthenogenetic reproductive events have caught 
the attention of evolutionary and conservation biologists 
since the absence of genetic recombination accelerates the 
accumulation of deleterious mutations in parthenogenet¬ 
ic individuals, which may also have implications in the 
management and conservation of the species (Ennos et al. 
2005, FIedrick 2007, Lampert 2008). In fact, partheno¬ 
genesis has been related to an increased incidence of abor¬ 
tive events, stillbirths, and embryonic and developmental 
abnormalities (e.g. anophthalmia, microphthalmia, en- 
cephalocoele and head foreshortening) in several species 
(Billy 1986, Booth & Schuett 2016, Allen et al. 2018). 

Obligate parthenogenesis is a biological process in 
which all individuals within a species reproduce asexually 
(van der Kooi & Schwander 2015, Booth & Schuett 

2016) . Among vertebrates, this reproductive strategy is re¬ 
stricted to Squamate reptiles, being reported in various liz¬ 
ard species and one snake ( Indotyphlops braminus ) (Wynn 
et al. 1991, Kearney et al. 2009, Booth & Schuett 2016). 
The occasional occurrence of parthenogenesis in individu¬ 
als of a species that normally reproduces sexually (i.e. fac¬ 
ultative parthenogenesis - FP) was first mentioned in birds 
in the late 1800s (Oellacher 1872). The phenomenon of 


FP has been reported in various species of major vertebrate 
groups including reptiles, birds and elasmobranchs (sharks 
and rays) (Booth & Schuett 2016, Harmon et al. 2016, 
Dudgeon et al. 2017, Ramachandran & McDaniel, 
2018). Most FP events were documented in captive females 
after long periods without contact with male conspecifics 
during their reproductive lifetime (Lampert 2008, Booth 
& Schuett 2016). However, parthenogenesis has recently 
been reported in both wild snake populations (Booth et al. 
2012, Calvete et al. 2018) and females housed with males 
(Booth et al. 2014), suggesting that its occurrence may 
be more frequent in vertebrates than previously thought. 
In addition, the reproductive viability of parthenogenetic 
offspring has been observed in some species, which high¬ 
lights the ecological and evolutionary significance of this 
reproductive strategy (Calvete et al. 2018). Nevertheless, 
the biological basis and mechanisms underlying partheno¬ 
genesis remain mostly unknown (Booth & Schuett 2016, 
Dudgeon et al. 2017, Shibata et al. 2017). 

Accurate identification and characterization of parthe¬ 
nogenesis in captive individuals of non-model species may 
provide important data towards understanding the fre¬ 
quency, causes, consequences and biological mechanisms 
of asexual reproduction among vertebrates (Fujita & Mo¬ 
ritz 2009, Neaves & Baumann 2011). The occurrence of 
FP has been described in at least six snake families, namely 
Boidae, Pythonidae, Viperidae, Acrochordidae, Colubri- 
dae and Elapidae (Booth & Schuett 2016, Shibata et al. 
2017, Allen et al. 2018). Among boid snakes (Boidae fam¬ 
ily), parthenogenesis has been recorded and validated by 
genetic analysis in Boa constrictor and Eunectes murinus, 
as well as in two species of the genus Epicrates (E. maurus 


© 2020 Deutsche Gesellschaft fur Herpetologie und Terrarienkunde e.V. (DGHT), Mannheim, Germany 
Open access at http://www.salamandra-journal.com 


83 




Correspondence 


and E. cenchria), closely related to the genus Chilabothrus 
(Booth et at. 2011a, 2011b, Kinney et al. 2013, Shibata et 
al. 2017). Various FP events confirmed by genetic analysis 
have also been reported in Pythonidae ( Python bivittatus, 
P. regius and Malayopython reticulatus, Groot et al. 2003, 
Booth et al. 2014), Viperidae ( Agkistrodon contortrix, Jor¬ 
dan et al. 2015, Calvete et al. 2018), Colubridae ( Tham- 
nophis marcianus and T. couchii, Germano & Smith 2010, 
Reynolds et al. 2012), and Elapidae ( Oxyuranus scutellatus 
and Acanthophis antarcticus, Allen et al. 2018). The appli¬ 
cation of molecular techniques greatly improved the char¬ 
acterization of FP events in wild and captive snakes. How¬ 
ever, to understand the evolutionary significance of this al¬ 
ternative reproductive strategy it is essential to increase the 
effort to detect evidence of new potential FP occurrences. 

Here, we report two interesting reproductive events re¬ 
corded in a captive female Cuban boa ( Chilabothrus angu- 
lifer) isolated from males for 11 years. These occurrences 
could be explained by two hypotheses: (i) long-term sperm 
storage from the last mating event or (ii) parthenogenetic 
reproduction. We applied molecular and histopathological 
methodologies to evaluate both hypotheses, providing the 
first evidence of facultative parthenogenesis in a Cuban boa. 

On September 20,2017, a 20-year-old captive female Cu¬ 
ban boa from the Zoo da Maia (Maia, Portugal) collection 
gave birth to a stillborn and deposited multiple non-em- 
bryonated eggs. This female, added to the zoological col¬ 
lection in 1999, has had no contact with a male since 2006 
when its mate died. Previously, in 2013, this female had de¬ 
livered a yellowish mass of non-developing eggs and a dead 
neonate, that has been preserved in 10% buffered formalin. 
The offspring from these two deliveries, one fresh and the 
other formalin-fixed, were analysed in the Histology and 
Anatomical Pathology Taboratory of Tras-os-Montes e Alto 
Douro University (UTAD). Tissues samples were processed 
for histopathology applying the routine staining techniques, 
with haematoxylin and eosin (HE), for light microscopy. 

The DNA isolation from the formalin-fixed specimen 
(neonate 2013) was carried out using the Quick-DNA 
Miniprep Plus Kit (Zymo Research), according to the man¬ 
ufacturer’s recommended protocol, with some additional 
steps before sample digestion. Briefly, a mixture of differ¬ 
ent tissues (liver, lung, gut and skin) was sliced into small 
pieces with a scalpel. Then the tissues were washed with 
Phosphate-buffered saline (PBS) during 24 h (the buffer 
was replaced twice). The DNA extraction from muscle tis¬ 
sues of the neonate born in 2017 was performed using the 
NZY Tissue gDNA Isolation kit (Nzytech). The mother’s 
DNA was isolated from her blood using the NZY Blood 
gDNA Isolation kit (Nzytech). Both extractions were per¬ 
formed following the standard protocols recommended by 
the manufacturer. 

Thirteen microsatellite markers previously character¬ 
ised for boid species were analysed: psat 1, psat 10, psat 
13, psat 24, psat 32, psat 36, Q25, Q34, Q35, Q36, Q37, 
55HDZ554, and 55HDZ617 (Ramanana et al. 2009, Tzika 
et al. 2009, Booth et al. 2011b, Reynolds et al. 2014). Pre¬ 
screening of microsatellite variations between mother and 


offspring samples was performed using high-resolution 
melting (HRM) analysis (Mader et al. 2008, Lombal et 
al. 2015). PCR amplification and melting acquisition were 
carried out using a QuantStudio 3 Real-Time PCR Sys¬ 
tem (Applied Biosystems). The reaction mixture was pre¬ 
pared to a 20 pi final volume containing 10 pi of MeltDoc- 
tor HRM Master Mix (Applied Biosystems), 5 pmol of each 
primer and 5 ng of genomic DNA. All PCR reactions were 
performed in duplicate. The amplification protocol was run 
as follows: 1 cycle of 95°C for 10 min; 40 cycles of 95°C for 
15 s, 6o°C for 1 min (fluorescence signal was captured at the 
end of each cycle); 1 cycle of 95 °C for 15 s, 6o°C for 1 min 
and then sequential temperature increments of o.o25°C/s 
with temperature ranging from 6o°C to 95°C, with contin¬ 
uous fluorescence measurements. The melting curve data 
were analysed using the QuantStudio Design & Analysis 
software v.1.4 (Applied Biosystems) and the High Resolu¬ 
tion Melting (HRM) Software v.3.0.1 (Thermo Fisher Sci¬ 
entific). The differences in melting-curve shapes allowed 
the characterisation of the microsatellite allelic variabili¬ 
ty. Forward primers of the microsatellites with variations 
among samples were labelled with 6-FAM to determine the 
genotypes using capillary electrophoresis. PCR amplifica¬ 
tions were performed to a total volume of 20 pi containing 
10 pi of 2 x MyTaq HS Mix (Bioline), 5 pmol of each primer 
and 5 ng of DNA. PCR thermal conditions were as follows: 
initial denaturation at 95°C for 5 min, followed by 35 cycles 
of 95°C for 30 s, 6o°C for 1 min, 72°C for 30 s and a final 
extension at 6o°C for 10 min. Amplified fragments were 
submitted to electrophoresis on an ABI PRISM 3130x1 Ge¬ 
netic Analyzer (Applied Biosystems) using the GeneScan 
500 LIZ size standard. Allele sizes were determined using 
Peak Scanner v.3.0.2 (Thermo Fisher Cloud). 

Our results support the first evidence of FP in the Cu¬ 
ban boa ( Chilabothrus angulifer). Macroscopically, both 
reptiles corresponded to fully developed stillborn snakes 
that showed no morphological alterations except in the 
2017 specimen (#2), which presented bilateral anophthal¬ 
mia (Fig. 1). Microscopic examination of organs showed no 
alterations and the inspection of the reproductive system 
confirmed both to be female. Of the 13 microsatellite loci 
screened, two markers (psat 13 and Q34) did not amplify 
or they generated non-specific PCR products. The high- 
resolution melting (HRM) analysis allowed the identifica¬ 
tion of four microsatellite loci (psat 10, psat 24, Q36 and 
Q37) with allelic variability between mother and offspring 
samples. No evidence of allelic variability was detected in 
the remaining loci for the samples analysed. These results 
were validated using capillary electrophoresis to determine 
allele sizes for all polymorphic markers and for two non- 
polymorphic loci (Table 1). Maternal heterozygosity was 
observed for polymorphic loci and a homozygosity pattern 
was obtained for the non-polymorphic loci analysed (Ta¬ 
ble 1). The offspring were homozygous for all microsatellite 
loci, always carrying an allele present in the mother (Ta¬ 
ble 1). The locus Q37 is a potential null allele in the neonate 
of 2013 since the amplification failed using different DNA 
samples and PCR conditions. 
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Table 1. Genotypes of mother and offspring characterised for the potentially parthenogenic Cuban boa (Chilabothrus angulifer) 


Individual 

psat 10 

psat 24 

Ci35 

Ci36 

Ci37 

55HDZ617 

Mother 

401/409 

215/217 

318/318 

220/228 

278/298 

185/185 

Neonate 2013 

409/409 

217/217 

318/318 

220/220 

- 

185/185 

Neonate 2017 

409/409 

215/215 

318/318 

228/228 

298/298 

185/185 


The allele homozygosity in the offspring, carrying only 
maternal alleles, was similar to previous FP cases reported 
for boid species (Booth et at. 2011a, 2011b, Kinney et al. 
2013, Shibata et al. 2017). In the Boidae family (for gen¬ 
era Boa and Epicrates), parthenogenesis was initially de¬ 
tected in all-female litters (Booth et at. 2011a, 2011b). Sim¬ 
ilarly, the two stillborn analysed in this study were both 
females. Accidental FP in captive individuals generally oc¬ 
curs after long periods of isolation from mates (Booth & 
Schuett 2016). The female in this study has had no con¬ 
tact with any conspecific male for 13 years now. Prolonged 
sperm storage has been documented in various snake spe¬ 
cies, but the longest time period of suspected sperm stor¬ 
age reported for a snake was seven years and six months 
(Magnusson 1979, Birkhead & M 0 ller 1993, Booth & 
Schuett 2011). In our case we described the lack of male 
genetic contribution to the offspring using molecular anal¬ 
ysis, thus excluding prolonged sperm storage. The high lev¬ 
els of homozygosity detected in the offspring are a charac¬ 
teristic of the parthenogenetic mode explained by terminal 
fusion automixis (Booth et at. 2011a, 2011b), as recently 
inferred for long-term captive copperheads (Agkistrodon 
contortrix ) and cottonmouths (A. piscivorus ) (Shibata et 
al. 2017, Allen et al. 2018). 

Genome-wide homozygosity of parthenogenetic off¬ 
spring maybe related to the development of malformations 
(Booth et al. 2012). Embryos and stillborn offspring with 
developmental abnormalities have been associated with 
parthenogenetic events in reptiles (Billy 1986, Booth & 
Schuett 2016, Allen et al. 2018). The morphological eval¬ 
uation of the Cuban boa neonate born in 2017 evidenced 
bilateral anophthalmia, a malformation found in parthe¬ 
nogenetic offspring of other reptile species (Billy 1986). 


In conclusion, we identified and described a probable 
first record of FP in the Cuban boa supported by speci¬ 
men history and molecular markers. This may have impor¬ 
tant ecological and evolutionary implications, including 
the improved understanding of this reproductive strategy 
and its frequency in captivity, and possibly in the wild, as 
recorded for some species (Booth et al. 2012, Fields et 
al. 2015, Calvete et al. 2018). The increasing number of 
reports on parthenogenic births in a wide range of snakes 
and other vertebrates also points to the evolutionary sig¬ 
nificance of this reproductive phenomenon, as yet poorly 
understood, making it a research field with high potential 
(Lampert 2008, Avise 2015). Therefore, we encourage zoo 
workers, veterinarians, curators, wildlife managers and re¬ 
searchers to pay attention and to investigate further evi¬ 
dence towards abnormal births in this species and related 
taxa, since these events may be easily overlooked. 
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